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Study of Cryogenic Propellant Systems for Loading 
the Space Shuttle* 


R. O. Voth, W. G. Steward, andW. J. Hall 

Computer programs have been written to model the liquid 
oxygen loading system for the space shuttle. The programs allow 
selection of input data through graphic displays which schematically 
depict the part of the system being modeled . The computed output 
is also displayed in the form of graphs and printed messages . Any 
one of six computation options may be selected. The first four of 
these pertain to thermal stresses, pressure surges, cooldown times, 
flow rates and pressures during cooldown. Options five and six 
deal with possible water hammer effects due to closing of valves, 
steady flow and transient response to changes in operating conditions 
after cooldown. 

Procedures are given for operation of the graphic display 
unit and minicomputer. 

Key words: Computer modeling; cooldown; cryogenic flow; stresses; 

surges; transient flow; water hammer. 

1 . Introduction 

This is a final report under NASA Contract No. CC-03420A. 


Loading of the space shuttle with liquid oxygen, using the existing 
ground support equipment at the Kennedy Space Center (KSC) , will require 
alterations in the current cooldown and loading procedures . Also minor modi- 
fications of the existing loading system may be required. Since acceptable 
loading rates, minimum loading times, and liquid capacities for the shuttle 
have not been fixed, the loading procedure and loading sequence has not been 
determined. Even after the shuttle capacities, etc. , have been determined, 

* This work was conducted at the National Bureau of Standards under the 
sponsorship of Kennedy Space Center - NASA . 


changes to the loading procedure will probably occur as experience is gained . 

Any changes in the cooldown procedure can cause excessive cooldown pressure 
spikes and/or excessive temperature induced stresses in components of the 

\ 

ground support system during the cooldown period. Also, flow rates of both 
liquid oxygen and vented gaseous oxygen during loading operations that follow 
cooldown are important. The study supplies the necessary mathematical and 
system analysis and hardware and computer software so that any proposed 
cooldown and loading sequence can be simulated; thus, the consequences 
and/or benefits of altering established procedures can be determined before 
a wet test of the system is initiated. The present study is for the liquid oxygen 
system . 

1 . 1 Objectives 

The objective of the study is to supply the KSC with the necessary 
analysis , hardware , and computer software to provide a computerized model 
of the liquid oxygen system. By means of these programs, cooldown and load- 
ing procedures can be checked mathematically to determine thermal stresses, 
cooldown pressure surges, cooldown times, and loading rates before a wet 
test of the actual system is conducted. This work is an extension of a recently 
completed study in which a mathematical model of the 35. 6 cm (14-inch) cross 
country piping system was developed for the calculation of cooldown pressure 
surges. Specific objectives are: 

( 1) Purchase of a graphical display unit with the hardware necessary 
to interface the display unit with an existing computer at KSC . 
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(2) Develop a mathematical model and computer program to 
calculate thermal stresses in major components of the ground support 
equipment during the cooldown period. The major components include 
the liquid oxygen pumps and other massive components such as pipe 
flanges. 

(3) Develop a mathematical model and computer program to calculate 
cooldown times for the various configurations of the ground support equipment. 

(4) Develop a mathematical model and computer program to 
calculate flow rates and pressures at selected points on the ground 
support equipment after the cooldown period. The flows and pressures 
shall be related to pump speed, supply tank liquid level, system con- 
figuration, and valve characteristics. 
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2. Graphic Display Computer Terminal 

2. 1 Description of Equipment 

The system consists of a graphical display terminal interfaced 
with the KSC time share computer through a 8000 word storage mini- 
computer and MODEM. Figure 1 shows the graphical display unit and 
minicomputer. 

Software for communication between the display unit and the KSC 
computer, routines for schematic displays, input-output, and graphing 
are programmed into the minicomputer. These routines are also stored 
on punched paper tape from which they may be reloaded if necessary. 

The programs and subroutines used in the calculations are catalogued 
in the KSC computer and are called as needed through the minicomputer 
software routines. 

2. 2 Description of the Minicomputer Software 

The programs contained in the minicomputer control the informa- 
tion flow between the graphical display unit and the KSC timesharing 
computer. Besides the control function, the minicomputer has the input 
schematics with the associated input data and programs to interpret the 
responses received from the KSC timesharing computer (during execution 
of the oxygen ground support computer program) in its memory. 

The communication between the minicomputer system and the 
KSC computer is carried over a high speed telephone data line. The 
MODEM properly conditions the signals for the high speed data 
transmission. 

Data is transmitted in various sized blocks. Each character of 
data in the block is accompanied by a parity bit and each block of data 
is accompanied by a parity character to make the parity of each word 
and parity of the entire block being transmitted correct. The receiving 
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station checks the parity of the incoming block of data and if the parity 
is correct, the transmission is accepted and an acknowledgment is 
returned to the sending station. If the parity of the received information 
is not correct, due to noise or other problems with the telephone line, 
the receiving station requests a retransmission. The sending station 
will then automatically retransmit the data. ; In some rareginstances, 
the retransmission request is lost and the data must be manually 
retransmitted. .'$'//'■ ' •' 

The minicomputer .software has two basic modes of operatibn 
using the graphical display keyboard. The first mode, termed the 
manual mode, allows conversation with the KSC computer in a manner 
similar to a normal teletype; except that the communication is at high 


“• speed. The second mode of operation, termed the automatic mode, is 
used to run the computer programs simulating the oxygen ground support 
equipment. Much of the communication between the minicomputer and 
the timesharing computer in the second mode is automatic, with no 
operator input required. ,• . .... 

During operation in either mode, the graphical display terminal 
will emit an audible "beep" whenever an acknowledgment is received 
-from the timesharing computer. This "beep" will inform the operator 
that a data block has been successfully sent and new information can be 
entered or a return can be expected from the timesharing computer. 

Until a successful transmission is affected, the keyboard of the graphical 
display unit is locked except to a control B , as explained in a later 
paragraph. 

Operation in the manual mode allows for normal log-in procedures 
and the use of any subsystems on the KSC computer. Fortran, basic 
or the editing subsystem, may be employed in the manual mode. 
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Since the graphical display -mini computer combination is connected 
to the Datanet-760 terminal of the KSC timesharing system, operation in 
the manual mode is very similar to the general operation of the Datanet 
system as described in the Honeywell manual [1971]. * Log-on procedures 
as described for the Datanet are also valid for the graphic display unit. 

During operation of the graphics unit, data entered on the screen 
of the graphics display are also entered into the memory of the mini- 
computer and when a carriage return is depressed, the data from the 
computer memory is transmitted to the timesharing computer. More 
than one line of data may be transmitted in a block if each line entered 
on the screen is terminated by a line feed and the entire block is termi- 
nated by a carriage return for transmission. 

Typing errors can be corrected by using a control A which causes 
the cursor to be backspaced once and lowered one line. The character 
typed in the lower position replaces the original character in the mini- 
computer memory and will be transmitted when the carriage return is 
depressed. Since the graphics display unit is connected to the Datanet- 
760 terminal, other characteristics unique to the Datanet system effect 
the operation of the display unit. As an example, the timesharing com- 
puter counts the lines of output during its data transmission and after 
sufficient lines have been sent to fill the screen of a Datanet-760 system, 
data transmission is terminated. The stop in data transmission is 
indicated by the ASCII code characters A * on the screen of the graphics 
unit. To receive additional data, the operator must transmit a blank 
line by depressing the carriage return on the keyboard. 

The automatic mode of operation interacts with the specific oxygen 
system simulation programs available in the Fortran subsystem of the 
timesharing computer. During automatic operation, the minicomputer 

* Identification of a manufacturer and a manufacturer's product has 
been necessary to make the results of this work meaningful and in no 
way implies a recommendation or endorsement by the National Bureau 
of Standards. 
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will request a run of a specified program, send the required data when 
an equal sign is received, and interpret and display returning data from 
the timesharing computer on the graphical display unit. Usually, the 
returning data will be displayed in graphical form. Since the oxygen 
simulation programs are in the Fortran subsystem, the KSC computer 
must be in the Fortran subsystem before the automatic mode of opera- 
tion is entered with the graphical display unit. The Fortran subsystem 
of the timesharing computer is indicated by an asterisk displayed on 
the screen. 

Several control characters are used during normal operation of 
the graphic display unit to help the operator during the execution of the 
various programs. A control character is used during communication 
between the graphic display unit and the minicomputer and need not 
be followed by a carriage return to effect a transmission. A control 
character is executed by depressing and holding down the control key 
on the keyboard and then depressing the desired alpha-numeric key. 
Thus, a control D means that the control key was depressed and held, 
and then the letter D was depressed. The following is a description of 
the control characters available for use on the graphic display unit. 

Control A: Control A allows the operator to correct a typing 
error before a block of data is sent to the timesharing computer. Any 
time before the carriage return is depressed, a change can be made 
to previously typed material. This change is made by depressing the 
control A an appropriate number of times until the cursor is located 
below the character to be changed, entering the new character and 
then proceding to enter subsequent characters in the line. 
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Control B: Depressing control B will cause the minicomputer to 

issue a $*$ BRK command to the timeshare computer and return the 
minicomputer to the manual mode of operation. The function of the 
$*$BRK command is described in the Honeywell Information Manual [1971] 
Control B can be used any time even during a transmission when the 
keyboard is locked to other input. 

Control G: Control G causes the minicomputer to enter the 
automatic mode of operation. Entering the automatic mode of operation 
should not be attempted unless the KSC timeshare computer is operating 
in the Fortran subsystem and an asterisk (*) is displayed as the last 
received character on the graphic display unit. 

Control P : Depressing the control P causes the minicomputer 
software to interpret data returning from the timesharing computer into 
graphical output. This option was supplied with the minicomputer to 
allow graphical output from any Fortran programs written in the manual 
mode. In order to graph the results of any program written in the 
manual mode, the output from the program must be supplied to the 
ASCII coded Fortran subroutine called PLOT currently available on the 
timesharing computer. Along with the data, labeling information can 
also be supplied to PLOT. The program DEMO currently loaded in the 
timesharing computer can be listed to obtain instructions on how to use 
subroutine PLOT. After writing a program for which a graphical output 
is expected the program is run in the manual mode, except that after 
all input data has been supplied and an acknowledge received as 
indicated by a "beep;" control P should be depressed. A graph of the 
returning data will then appear on the graphical display unit. 
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2. 3 Instructions for Log-In 

The log-in routine for the minicomputer is accomplished in the 
manual mode. The procedure is to turn on the Modem (switch located 
inside the front cover), the minicomputer (key switch on front panel), 
and the graphics display unit (switch located on the front of the pedestal 
below the keyboard). The line -local switch on the graphics display 
keyboard should be set to line. After a short warm-up period, the 
screen on the graphics display unit will light up. The page key on the 
display unit should be depressed to prepare the screen for reception 
of messages. The starting address of the software program in the 
minicomputer is 00002 in octal (switch 14 on the top row of toggle 
switches up, all others down). After setting the top switches to the 
starting address, depress RESET and then START on the momentary 
contact control switches from the lower row of switches of the mini- 
computer. The timesharing computer can now be dialed on the tele- 
phone. After the computer carrier signal (a continuous tone) is 
confirmed, the telephone data switch can be actuated. As soon as the 
minicomputer software confirms the computer interconnection, the 
words PLEASE LOG IN will appear on the graphics display screen. If 
these words do not appear, a fault has developed in the minicomputer 
software, or the telephone interconnection was not achieved. The 
program should be reloaded as described in appendix C. The normal 
log-in procedure for the Datanet 760 system as described in the 
Honeywell manual [1971] can now proceed. A carriage return should be 
supplied after each input to transmit the information to the timesharing 
computer. Figure 2 shows the appropriate log-in sequence as it would 
appear on the screen of the graphics display unit. Computer supplied 
messages are underlined. 
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PLEASE LOG IN 

$*$MD PASSWD, 26, DAC.TSS 

CONNECT REQUEST SENT TO CP 
KSC 635 AT CIF ON 03/06/74 AT 10. 281 CHANNEL 4424 


USER ID - 
3MX002 

PASSWORD— 

NBS 

9 BLOCKS FILE SPACE AVAILABLE 
**TSS SIZE 54K CARDIN 25K 


SYSTEM ? 
FORT 

OLD OR NEW- 
NEW 

READY 

* 


Figure 2. Log-in sequence. 
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2. 4 Instructions for Initiation of Computation in the Automatic 
Mode 

In order to perform calculations in the automatic mode using the 
oxygen simulations programs, the Fortran subsystem of the time- 
sharing computer must be available. On figure 2, the Fortran sub- 
system was entered by entering Fortran for the required system and 
new for the type of program. After an asterisk (*) appears on the 
screen as the last reception, control G can be depressed to enter the 
automatic mode of calculation. After initiation of the automatic mode, 
the universal figure (fig. 3)* will appear on the screen of the graphics 
terminal. This figure is used to define the general characteristics of 
the piping system. 

Pipe lengths under Entry 1 in the General Input Data may include 
equivalent length to account for the added restrictions of pipe fittings, 
filters, or other restrictions. If this is done, the pipe wall thickness 
entered into the specific input data of option 3 must be altered corre- 
spondingly in order to provide the correct mass of material to be 
cooled down. Entry 7, "Supply Ullage Pressure" requires some clari- 
fication. For computation options, except option 6, there is no differen- 
tiation between the ullage pressure during storage and during the trans- 
fer process. For option 6, the pressure specified in the general input 
data is taken to be the ullage pressure during storage and it is assumed 
that the supply liquid becomes saturated at that pressure. The specific 
input data for option 6 allows a "Supply Ullage Pressure During Transfer" 
to be selected which may be different from the ullage pressure during 
storage. 

When the figure appears, the input variables will contain the 
numbers used in the previous run of the program. Any or all of the 


The figures shown here are actual photographs of the graphic display 
as they will appear during operation of the system. 
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input numbers may be changed. A change is accomplished by moving 
the horizontal cursor to the line containing the number to be changed, 
pressing the "RUB OUT" key, then typing the new number in the blank 
space which will appear. After pressing "RETURN" the operator may 
move the cursor and make further changes. Even though the schematic 
will not change, the computer programs will accept zero for lengths 4 
and 5 (L4, L5) to eliminate the vertical rise if desired and will accept 
a single pipe length equal to the total pipe length and a single diameter. 
At the upper right on the figure, the current date and time will appear. 
This date and time will remain the same and will appear on all subse- 
quent figures to be used as an identification number to correlate hard- 
copy output. After all input parameters have been adjusted to acceptable 
values, and the ALT MODE key is depressed, the system will proceed 
to the next page. 

This page (fig. 4) allows a choice of the computer program to be 
run with the previous input data. After this choice is made by placing 
the option number in the appropriate blank and depressing ALT MODE, 
the terminal will display the first page of that option. The first page 
of each option, plus as many as two more pages will appear, allowing 
for more specialized input data. Typical display figures are shown in 
the following sections. After the input data has been adjusted as 
previously described and ALT MODE depressed for each input figure, 
the output of each option will appear after a short period of time. 

During this time interval, several "beeps" will be heard as input data 
is supplied to the timesharing computer. In some instances, when 
multiple output is expected, a choice of output curves will be made 
available. Each output curve will fill the entire screen. When the 
output curve appears, it will be possible to make a hard copy and/or 
digitize any point on the curve by placing the vertical cursor on the 
desired point and depressing the space bar. 
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Restarting with the input display is accomplished by depressing 
ALT MODE. At this time the operator may choose another computation 
option with the previous input values, he may change the input, or he 
may terminate computation by selecting option 7 from the second page 
(RETURN TO MONITOR). Of course, if the operator chooses to start 
calculations on a new system, a new date and time will appear on all 
figures, allowing him to maintain correlation between the input and 
the output figures. 

2.5 Reloading the Minicomputer Programs 

Several paper tapes have been supplied with the system. Some 
of the tapes were supplied by the manufacturer of the minicomputer 
and can be used to debug the hardwired functions of the minicomputer. 
The function of the manufacturer's supplied' tapes are described in 
accompanying literature. Other paper tapes can be used to reload the 
existing minicomputer program. ’ The paper tape loading procedure, 
using an appropriate teletype unit, is described in Appendix C. 
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Figure 4. Selection of computation options. 



3 . Computation Options — Cooldown Period 

3.1 Option 1 - Thermal Stresses 

Calculates stresses due to rapid cooldown of flanges, or bowing of 
pipeline due to stratification in the case of low flow rates . 

3.1.1 Applications 

This computation determines low and high limits of flow rates during 
cooldown from the standpoint of the resulting stresses in the piping system 
and flanged joints. For very high flow rates and rapid cooling excessive 
thermal stresses may be created in flanges due to the large radial tempera- 
ture gradients. On the other hand very low flow rates might result in 
stratified flow with liquid lying along the bottom of the pipe . The resulting 
circumferential temperature gradients could cause a pipeline to bow significantly . 

3.1.2 Input Data Entries 

Section 2 . 3 describes the general log-in procedure for the graphic 
display system. After the initial log-in has been accomplished, and Option 1 
has been selected, a display similar to figure 5 appears. The input required 
for this option in addition to the general input data of figure 3 is shown in 
figure 5 . These are the outside diameter of the pipe flange, the outside and 
inside diameters of the pipeline and the material of the pipeline . 

3 . 1 . 3 Output 

The output for this calculation is alpha numeric only but will be dis- 
played on the screen just as other output data. The values given 
will be the maximum cooldown flow rate which will allow the stress in the 
flanges to remain below the yield point, and the minimum flow rate which will 
prevent stratified flow from occurring with the given input conditions . 

17 
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Figure 5. Option 1 - Thermal stresses input display. 



3.2 Option 2 - Valve Opening Surges 

Surges due to rapid vaporization when a valve is opened preceding a 
free wheeling pump, or water hammer as the liquid impinges on a free wheeling 
pump . 

3.2.1 Applications and Limitations 

This computation applies to the opening of a valve preceding a warm 
section of pipe and a pump . It computes pressure surges which might develop 
from one of two causes: (1) rapid vaporization in the warm pipe (model 

discussed in Appendix B) , or (2) a water hammer effect due to impingement 
of liquid on a free wheeling pump. When the vaporization surge 
occurs, it is accompanied by a slowing or reversal of the liquid flow so that 
the water hammer would not occur or would be greatly diminished. The 
program estimates the time at flow reversal due to vaporization and the time 
of liquid impingement at the pump in order to determine which of the two kinds 
of surges to compute . Details of the computation of these pressure surges and 
a computer listing of program VOS are given in Supplement A (1974), The 

derivation of the pressure surge correlation is given in Appendix B. 

3. 2. 2 Input Data Entries 

Section 2.3 describes the general log-in procedure for the graphic 
display system. After the initial log-in has been accomplished and the 
Option 2 has been selected, first a display similar to figure 6 appears 
and after ALT mode is depressed figure 7 will appear. If desired, input 
variables may be examined and changed in the display as explained in 
Section 2. 4. The variables in figure 6 and 7 apply only to Option 2 and 
are in addition to those of the general input data of figure 3. 
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Figure 6. Option 2 - Valve opening surges input display. 
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Figure 7. Option 2 - Valve opening surges input display 
(continued). 



3. 2. 3 Output 

The output of the computation is a graph of peak surge pressure 
as a function of valve opening time. A typical example is shown in 
figure 8. As the valve opening time is shortened, the peaks increase 
to a maximum amplitude at a critical value of the opening time, as 
discussed in Appendix B. Further decrease in the opening time has no 
effect on the peaks. The maximum pressure and corresponding critical 

opening time are indicated in the graph. If the water hammer surge is 

/ 

applicable, no graph is displayed, but a single amplitude of the pressure 
spike is given. 
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Figure 8. Option 2 output. 



3. 3 Option 3 - Cooldown Time - length cooled down, flow rate, 
and pressure drop across the pump versus time. 

3. 3. 1 Application 

As a cryogenic liquid flows or is pumped from a reservoir into 
a warm pipeline, vapor forms at the leading end of the liquid stream. 
Frictional resistance to the flow of this vapor is the principle limiting 
factor in the rate of cooling pipelines of large length- to- diameter ratios. 
The slow moving liquid is evaporated and vapor is warmed to near the 
initial temperature of the pipe in a short distance relative to the length 
of pipe. When the pipeline is uninsulated, the heat transfer to the cold 
part of the pipe adds to the vaporization rate and thus further retards 
the progress of the liquid. The cooldown model based on this cold 
front' concept and listing of the computer subroutine (CDT) are discussed 
in Supplement A (1974). The heat-mass transfer correlations are 
discussed in Appendix A. 

3. 3. 2 Input Data Entries 

Section 2. 3 describes the log-in procedure for the graphic display 
system. After the initial log-in has been accomplished and option 3 
"Cooldown Time" has been selected, a display similar to Figure 9 
appears. If desired, input variables pertinent to the cooldown time 
computation may be viewed and changed on this display as explained 
in Section 2. 4. The variables in figure 9 apply only to Option 3 and 
are in addition to those of the general input data of figure 3. As 
discussed in Section 2. 4, if equivalent pipe lengths have been added to 
account for line restrictions, the pipe waill thickness should be decreased 
correspondingly such that the total mass of material remains unchainged. 
Concentrated masses such as valves, pumps, and fittings should be 
compensated for by specifying additional pipe wall thickness provided 
the calculated cooldown time desired is to include the time needed to 
completely cool these concentrated masses. 
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GENERAL INPUT DATA 65/29/73 14.464 
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Figure 9. Option 3 - Cooldown time input display. 



3.3.3 Output 

The outputs of the computation are: (1) Length Cooled Down versus 

Time, (2) Flow Rate of LOX versus Time, and (3) Pressure Drop across the 
Pump versus Time. Typical examples are shown in figures 10, 11, and 12. 
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3 . 4 Option 4 - Pump Start-up Cooldown Surges - surges due to rapid 
vaporization when a pump is started upstream from a warm line. 

3.4.1 Application 

The computation applies to the starting of a pump in a line which is 
partially or totally warm. By the correlation discussed in Appendix B it 
computes pressure surges due to rapid vaporization when liquid enters a warm 
line . 

3.4.2 Input Data Entries 

Section 2 . 3 describes the general log-in procedure for the graphic 
display system. After the initial log-in has been accomplished and the 
Option 4 has been selected, a display similar to figure 13 appears. If de- 
sired, the input variables may be examined and changed in the display as 
explained in Section 2 .4 . Only two variables , length of cold line and final 
pump RPM , are shown in this display. The remaining values for this cal- 
culation were taken from the general input data of figure 3. 
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3. 4. 3 Output 

The output of this computation is a graph of peak surge pressure 
as a function of pump start-up time. A typical example is shown in 
figure 14. As the pump start-up time (time to bring the pump from 
rest to operating speed) is shortened, the amplitude of the peaks 
increases until a maximum is reached at the critical value of the 
start-up time. Further decrease in the start-up time has no effect 
on the peak amplitude. The maximum pressure and corresponding 
critical pump start-up time are indicated in the output graph. 
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4. Computation Options, Post-Cooldown Period 

4. 1 Option 5 - Water Hammer - pressure spikes due to water 
hammer when valves are closed in a flowing liquid stream. 

4. 1. 1 Application 

This computation applies to the case in which a downstream valve 
is closed against a flowing liquid stream. Deceleration of the liquid 
due to impingement in the closing valve can lead to a sudden large 
pressure spike. The amplitude of the pressure spike is a function of 
the valve closing time. Details of the computation and computer 
listing are given in Supplement A (1974). 

4. 1. 2 Input Data Entries 

Section 2. 3 describes the general log-in procedure by which the 
Option 5 is selected. After this selection has been made, a display 

similar to figures 15 and 16 will appear in sequence. In the display, 

t: " 

if desired, the input variables may be examined and changed by the 
procedure explained in Section 2. 4. 
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Figure 15. Option 5 - Water hammer input display. 
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4. 1. 3 Output 

Typical graphs of the amplitude of water hammer pressure-<spikes 
as a function of the valve closing time is shown in figure 17. 
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4. 2 Option 6 - Cold Flow. Steady flow operation and transient 
response to changes in operating conditions. 

4. 2. 1 Application 

This computation applies after the transfer line has reached 
operating temperature. It deals with response to changes in operating 
conditions but does not deal with cooldown transients. The latter are 
treated separately in Section 3. The flow rate and liquid discharge 
temperature (or two phase quality) are calculated as functions of the 
pump speed, ambient air conditions, valve settings, and pipe geometry. 
The accumulated liquid, change in supply liquid level, and pump discharge 
pressure are also computed. Pump speed and valve settings are allowed 
to vary with time. Details of the computation and computer program 
listing are given in Supplement A (1974). 

4. 2. 2 Input Data Entries 

Section 2. 3 describes the general log-in procedure for the graphic 
display system. After the Option 5 is selected, a display similar to 
figure 18 and 19 will appear. If desired, input variables may be examined 
and changed in the display by the procedure explained in Section 2. 4 
except in the table of numbers showing time, valve opening, and pump 
speed, where the vertical and horizontal cursors must be positioned 
before depressing the space bar. 

Entry 1 allows a supply ullage pressure during transfer to be 
selected which may or may not be different from the ullage pressure 
during storage. (The storage pressure is taken as that pressure which 
was selected previously under General Input Data. ) Entry 1 along with 
entry 2 relates to the subcooling of the inlet liquid. This subcooling 
enters into the calculation of discharge liquid temperature, or fraction 
of vapor being discharged. Entry 4 allows for the selection of a schedule 
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of changes in operating conditions. That is, the pump RPM or the 
discharge valve setting may be specified at six times. RPM and 
valve setting change linearly from one specified time to the next. 

The program carries out the computation to ten seconds beyond the 
last entry in the schedule. 

The third entry under Cold Flow allows the selection of computed 
quantities to be displayed. The desired quantities are selected as 
numbers 1 through 6. Striking ALT MODE initiates the calculation, 
and the next display is the selected output graph. v 
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90 0.50 

100 V 0.50 

300 0.50 

Figure 18. Option 6 - Cold flow input display. 





4. 2. 3 Output 

Typical output graphs are shown in Figures 20 through 25. The 
following variables may be plotted as functions of time: 


Fig. 

20. 

Flow rate 

Fig. 

21. 

Total accumulated liquid discharged into the receiver 
after subtracting off the quantity vaporized 

Fig. 

22. 

Pump discharge pressure 

Fig. 

23. 

The line pressure ahead of the discharge valve 

Fig. 

24. 

Final vapor fraction 

Fig. 

25. 

Fluid discharge temperature . 
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STEADY OPERATION AND TRANSIENT RESPONSE 
LIQUID SUPPLY DEPLETED BEFORE FINAL TIME 
DISCHARGE TEMP AT A FLOW OF 10443 GPM WAS 162 R 
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Figure 20. Option 6 output. 





STEADY OPERATION AND TRANSIENT RESPONSE 
LIQUID SUPPLY DEPLETED BEFURE FINAL TINE 
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Figure 22. Option 6 output (continued). 
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STEADY OPERATION AND TRANSIENT RESPONSE 
LIQUID SUPPLY DEPLETED BEFORE FINAL TIME 




48 


Figure 24. Option 6 output (continued). 
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Figure 25. Option 6 output (continued). 
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Appendix A. Heat-Mass Transfer Correlations for Uninsulated 
LOX Transfer Lines 

Expressions for heat flux from the atmosphere to uninsulated 
LOX lines are needed for estimations of fluid enthalpy rise along the 
line. This enthalpy rise in turn determines whether the fluid being 
discharged is subcooled or is saturated and partly vaporized. A second 
application is in the computation of cooldown time. During cooldown 
the heat from the atmosphere adds to the vaporization rate due to 
cooling of the pipe walls and thus lengthens the cooldown time. 

After the pipeline has cooled down to a stable operating tempera- 
ture, the resistance to heat transfer is primarily in the outside film - 
that is, between the pipe wall and the ambient air. The heat flux then 

depends on the ambient air conditions and the build-up of frost on the 

2 2 

transfer line. A minimum heat flux of about 0. 076 W/cm (240 BTU/hr ft ) 
is expected in still dry air after the maximum thickness of insulating frost 

is established. Heat flux under the worst conditions should be less 
2 2 

than 0. 95 W/cm (3000 BTU/hr ft ). These values lie within the nucleate 
boiling range for liquid oxygen (Brentari, et al. [1965]) therefore, film boiling, 
with its large wall-to-LOX temperature differences, should not be 
necessary to consider. 

The heat flux after steady state is reached depends on the frost 
thickness and surface temperature. These factors are at least partially 
self-regulating in that the frost tends to build up to a stable thickness then 
begin to slough off. A frost surface temperature drop causes additional 
frost thickness to build up, increases the insulating effect, and thus 
re-establishes the higher frost surface temperature. This regulation 
is imperfect, however, and determination of this frost surface tem- 
perature is probably the most questionable factor in the heat flux pre- 
diction. The value of frost surface temperature calculated by Holten [1961] 
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and measured by Richards et al. [1962] for liquid nitrogen cooled 
surfaces -- corrected for the higher temperature of LOX -- was 
used in the present calculations and appears to give reasonable pre- 
dictions judging by comparisons with the liquid nitrogen data of 
Richards et al. as well as Ruccia and Mohr's [I960] LOX data. Our 
method over-predicts Holten's LOX data by 20 to 30%; however, the 
heat flux in those experiments may have been reduced from the free- 
air values by the partially sheltered placement of the test vessels. 
Holten's correlation fits his own data but underpredicts the other data. 

Heat flux occurring after cooldown, but before the frost has 
established a stable thickness, is considerably higher than the steady 
state values. This initial heat flux is predicted satisfactorily by 
assuming the surface temperature equal to the liquid saturation tem- 
perature. An empirical curve fitting was used to obtain the variation 
of heat flux with time between the initial and steady state values. 

Our developments of the correlations and comparisons with 
data are summarized in the following sections. 

Al. Mass Diffusion Heat Flux 

In heat transfer to surfaces at cryogenic temperatures con- 
densation of water vapor (or air) can along with convection and radia- 
tion, comprise a significant part of the total heat flux. Where con- 
densation is present with forced convection, it is assumed that the 
convective effect is not altered by the mass transfer and the heat 
fluxes resulting from the two effects may be added. Radiation is the 
third component to be added. 

The heat transfer effect of mass diffusion is due to the giving 
up of sensible heat and latent heats of vaporization and fusion by 
the water at the condensing surface. The mass diffusion heat flux is 
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(Q /A ) dif = w( I a .I.) 


(Al) 


w = rate of condensation per unit area 

I = enthalpy of the water vapor in the ambient air 

I. = enthalpy of the water in the condensed state 

I -I. as 2990 J/g (1287 BTU/lb). 
a i _ 

In terms of the diffusion coefficient 


w = h 


4 


wa 


dif \RT 


wi V 

RT i r 


wa 

l dif RT 


= h H p 
dif a 


(A2) 


dif 


wa 


wi 


a i / a 

= mass diffusion coefficient 

= partial pressure of water vapor in ambient air 

= partial pressure of water vapor at the condensing 
surface 


R = gas constant 

T = temperature ‘ 

(P,R,T in consistent units) 

H = specific humidity 

p = density of ambient air (due to the approximate 

- : nature of this calculation a standard value of 

0.0012 g/cm^ (0.075 lb/ft^) was assumed). 

: • The mass transfer of water vapor is treated through an 
empirical analogy between heat and mass transfer. The Chilton - 
Colburn analogy (Rohsenow and Choi [ 196l])gives 


dif _ ,2/3 

Tr~ ( Sc > 


p C V 
a pa 


(Pr) 


2/3 


or 


dif 


P c 

a p a 


(*■) 


2/3 


(A3) 
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where C = specific heat of air at constant pressure 

Pr = Prandtl number for air, C \i f K 

P 

= viscosity 

K = thermal conductivity 
Sc = Schmidt number = p. / pD 
, D = mass diffusivity 
V = velocity. 

From Al, A2, A3 


< Q/A > di f 





(A4) 


If it were not for mechanical failure of the frost, one would 
expect the outer surface to reach approximately 273 K (32° F); however, 
the frost surface continually sloughs off, exposing a colder layer. 
Therefore, as the water vapor diffuses toward the cold surface, it 
encounters a temperature at which it condenses before it actually 
reaches the surface, and much of the condensate falls or is swept 
away without attaching itself to the surface. Then some of the heat 
given up by the water vapor is indirectly transmitted to the pipe surface 
through the air convection currents. 

Holten [1961] calculated from his heat transfer data a steady 
state L.OX container frost surface temperature, T^, of about 210 K 
(- 80° F). Exactly the same figure can be derived from the frost surface 
temperature of 204 K measured by Richards, et al. (1962) for a liquid 
nitrogen test section if one assumes that the fractional temperature 
drop across the frost layer, is the same for LOX as for LN^ . That is; 


AT /AT = AT /AT 

frost, LOX' over-all, LOX frost, LN^ over-all, LN^ 

(T - 90)/(300 - 90) = (204 - 76)/(300 - 76) 

T = 210 K. 


The average film temperature is (210 + 300) /2 = 255 K (0°F). 
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The following properties are evaluated at 255 K (Rohsenow and 
Choi [1961]). 


for air Pr = 0. 722 and C = 1. 0 J/(g * K) (or 0. 239 -rr — ~ ) 

pa lb t 

for water Sc = 0.6. 

Substitution of these constant values gives 



= 3.38 (10) 3 hH (h [=] 


• :j W 

cm 2 K 


, , g (H 2 0) 

H f =1 — — ■ 

g (air) 


or 


(Q /A) 



0.868 hH'(h [=] 


BTU ; 
hr-ft 2 °F 


H[=] 


grain; 

lb 


). (A5) 


A2. Convection Heat Flux 

A2. 1 Natural Convection (still air) 

McAdams [1942] gives the following expression for natural 
convection across single horizontal tubes, 


Nu ' = h — = 0.47 [Gr Pr] 0,25 " ' ' - T (a6) 


Gr : is the Grashof number evaluated at the film temperature 


Gr = 


where 0 = thermal coefficient of volume expansion . <■. 
The group in parenthesis is given by Rohsenow and Choi 

for air at 0 °F. . ■ . ' 'V 


P 2 gR 

2 


= 267 


3 

cm K 


( 


or 4. 2( 10) 


3 

ft °F 


) 


- ~'p \ D 3 AT 
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From the same source 
Pr =0.722 


K 


= 2. 25(10) 


-4 W , BTU % 

cm K (OT 0-013 hr ft 0 F )- 


Substitution of the constant values into (A6) yields 

,1.25 


m 


c cm 


W „4 AT' 

— = hAT = 3.94(10) (AT [ = ] K, D[ = ] cm) (A7) 


D 


or 


(*•). 


BTU „ AT 
=- = 0. 255 


1.25 


c hr ! ft 
From (A5) we get 


D 


0. 25 


(AT [ = ] "F, D [ = ] ft).. 




0. 25 


H 


AT [ = ] K, D [ = ] cm, H[=] 


-1 £ 


(A8) 


or 


g, 

.0. 25 


H 


(2.) uSS* = 0 .221&f-) 

' A / di£ hr £t 2 \ D / 

AT[ = ] "F, D[ = J£t,H[ = ] 8^2! . 


A2. 2 Transition Velocity - Natural Convection to Forced 
Convection 

The transition velocity is obtained by equating Nusselt number 
for natural convection to Nusselt number for forced convection. For 
forced convection in the Reynolds number range 1,000 to 50,000 
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Nu = 0. 26 Re 0 ' 6 Pr°* 3 


(Me Adams [1942]) 


(A9) 


By equating (A6) and (A9) we get 


Re = 


pVD 


_ r £iii 
~ L 0 - 26 


47 _ 0. 25 Pr 
Gr 


0. 25 -i 1/0.6 


Pr 0 ' 3 J 


Substitution of the constant values listed previously gives 


or 



= 3.18 D°' 25 AT 0,417 (D [ = ] cm, AT [ = ] K) 

690 D 0,25 AT 0,417 (D [ = ] ft, AT [ = ] °F). 


For D = 30.5 cm (1 ft) and AT = 90 K (162 F) 


V — 49 or 1 . 76 Km/hr or 1. 1 mi/hr. 

s 

Therefore, if the wind velocity is less than approximately 
1 mi/hr, the natural convection expression will be used. If V is 
greater than 1 mi/hr, one of the following forced convection expres- 
sions will be used. 

A2.3 Forced Convection 

Equation (A9) fits experimental data given by McAdams up to 

5 

Re 2:10 . This corresponds to a wind velocity of about 12. 1 Km/hr 
(7.5 mi/hr). Equation (A9) diverges rapidly from the experimental 

5 

data above Re = 10 so that at 24. 1 Km/hr (15 mi/hr) (A9) underpre- 

5 

diets the data by 40%. For Re >10 (V > 7.5 mi/hr) we have developed 
the following equation, assuming the same Pr dependence: 
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Nu = 0.00332 Re 


(A10) 


0.982 0.3 

Pr 


A2. 3. 1 Heat Flux for 1. 76 < V < 12, 1 Km/hr 
(1 < V < 7. 5 mi/hr). 


From (A9) we get 


h = 0. 26 


K / pVD \ 0,6 
D { n ) 


(Pr) u ‘ J . 


Substitution of the constant values yields 

r 0. 6 


(-sr) 


-4 V 


cm 


1.75 (10) — q-j ( V [ = ] — ^ , D [ = ] cm) (All) 


D 


or 


( BTU \ * 8 

2 = °* 608 TT ( V [ = J ft / s * D [ = ] ft). 

hr ft °F / D 0,4 


From (All) 


r 0. 6 


(a") ~ V = 1 ; .754 (10f 4 AT (V[ = J-^, £) [=] cm 


c cm 


AT [ = ] K) 


or 


M.f, 


r 0. 6 


BTU = 0.666 ^r -7 AT (V[=] ft/s. D [=] ft. AT = ° F) (A12) 


ft D 

From (A5) and (All) 


0.4 




r 0. 6 


cm 


dif cm 


W V' 

T = - 59 °-(rr H (v AT{=]K, H[=]5- ,D[=]cm) (A13) 

^ a or 


D 


g 


or 




BTU 


r 0. 6 


dif hr ft 


= 0.5 78 H (V [ = ] ft/s, AT [ = ] X> F , H[=]^S£, 


lb 


D[ = ] ft) . 
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A2. 3. 2 Heat Flux for V > 7. 5 mi /hr 


From (A10) we get 


h = 0.00332 


K_ / pVD \ 

D W j 


0.982 


(Pr) 


0.3 


Substitution of the constant values gives 
W 


'em K / 


, v . 982 

75 <10) — ols" < v t = Jv-> D M cm > 


D‘ 


t /— ) = 8.428 (10) 
'hr ft K / 


-5 V 


0.982 


ft 


.0.018 


(V[ = J^. D[ = ]ft). 


Due to the small exponent of D, and the value of D being on the order 
of 30.5 cm (1 ft), 

0.018 


\ 3o * 5 / 


— 1.0. Then, 


(-?-)*■ 

\cm K / 


h( = 4.72 (10f 6 V 0 * 982 (V [ = ] — ) 

s 


(A14) 


or 




4. 72 (10)" 6 v 0 - 982 AT. 


From (A5) and (A14) 


(x) = 0.123 V°- ,82 H (V[=]2p-, H[=] £) 

' ' dif cm B 


(A15) 


( 


x) - BTlJ A = 0.2272 v 0, 982 H (V [=]— , H[ = ] ^v^ S ). 
A/ , 4 s 


'dif hr ft 


lb 
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A3. Radiant Heat Flux 


The heat flux due to radiation from the ground and sky is 


x = oF [< - T f )• 

For diffuse radiation from a distant enclosure the emissivity factor 
F may be taken as the emissivity of the frost surface. Holten [1961] 
and others assign values of abput 0.9 to this factor. The Stefan -Boltzman 
constant CT is 

- 1 2 W 

d = 5.67 (10) Y~T~ 

cm K 

or 

0. 174 (10) -8 BTU/hr ft 2 °R 4 . , 

T. = 210 K (378°R) as before. Substitution of these values give 
frost 

= 5.1 (10)- 12 (T 4 - 210 4 ) (T [ = ] K) 

• s’ „ / . 3» 



_Q 

A 


/ BTU \ 
\hr ft 2 °F / 


= 0.15 (10)' 


(T a " 


378 ) (T [ = ] °R) . 


In open sunlight the mean heat flux to an essentially black sur- 
face is approximately 


— = 0.068 Wcm 2 

A 

• p - ' 

= 215 BTU/hr ft 2 , 

A 

P 

assuming an atmospheric transmittance of 0.5. A^ is the projected 
area (seen by the sun). In terms of the pipe surface area A: 
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ttDL 

DL 


= tt . 


Therefore 


£. = °!°6£ = 0.0216 w/cm 2 or 68.4- — -- . 

A " hr a 2 

To avoid unnecessary complication, we assume 50% sunlight. 

the total radiant heat flux is 


i nen 


(H 


ad cm 


W -12 4 4 

— r- = 5.1 (10) (T - 210 )+ 0.0216 (T [ = ] K) (Al6) 

3> 3. 


(r) — 

' ' rad hr ft 


BT ^- = 0. 154 <10)" 8 (T 4 - 378 4 ) + 34. 2 (T [ = ]R). 

C, 3 > cL 


A4. Conversion of Relative Humidity to Specific Humidity 

Specific humidity H may be expressed in terms of the total 
(barometric) pressure B and the vapor pressure as 


H = 


gH 2 Q 


1. 61 (■#--!) * air 


or 


H = 


7000 


1.61 - D 

v 


grains 

lb 


The vapor pressure corresponding to the dry bulb temperature 

T °F is, by the Antoine equation (Lange [1961]) 
a 


= exp 


^18. 
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7254.4 
391 + T 


) 


mm Hg . 


(A17) 
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From the definition of relative humidity, 


P = H P, 
v r d 


Then, 


H = 


1 

1.61 (B/H r P d - 1) * 


An average value of B = 760 mm Hg yields 


H = 


1224 

H 

r d 


- 1.61 


£ 

g 


or 


(A18) 


grains 

lb 


(P from eq (A17)). 
a 

A5. Comparisons With Experimental Data 

Calculations of heat flux to uninsulated liquid and liquid oxygen 
transfer lines and vessels are compared with experimental data from 
the following sources: 

Author Fluid 

Richards et al LN 
(1962) 2 

Ruccia, Mohr LOX 

(1960) 

Hoi ten LOX 

( 1961 ) 


Apparatus 

Horizontal cylindrical test section 5.08 cm 
OD x 15. 2 cm long (2 in OD x 6 in long) in a 
wind tunnel. Controlled velocity and humidity. 

Vertical cylindrical tank 45. 7 cm OD x 6.71 m 
(18 in OD x 22 ft) high. Engine driven propeller 
provided wind. Humidity not controlled but 
measured. 

Spherical tanks 13. 7 cm OD ( 5-3/8 in). 

Natural convection. 


H = 


7000 

1224 

H 

r d 


1.61 
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A5.1 Comparison with Data of Richards et al. [1962] 

Since liquid nitrogen was used in these experiments the value 
of T was 204 K (368 °R) and AT was 95.6 K (172 °F) (Sec Al). The 
average frost thickness was approximately 0.952 cm (3/8 in), so that 
effectively 1.91 cm (3/4 in), or 1.91 cm (0.0625 ft) is added to the 
outside diameter. That is, D = 5.08 + 1.91 = 6.99 cm. Equation 
(A12) becomes 

(Q/A) = 0.00770 V 0,6 (V [ = ] cm/s) (A19) 

C u 

cm 

or 

(Q/A) j - = 206.6 V°- 6 (V [ = ] ft/s). 

C hr ft 

Equation (A13) becomes 

(Q/A) dif — ^2~ = - 272 V ‘ 6 H (V [ = ] cm/s , H [ = ] g/g) (A20) 
cm 

or ' - 1 ' 

(Q /A) d . f = 1. 042 V 0 * 6 H (V [ = ] ft/s, H [ =] grains /lb). 

And the heat flux due to radiation without sunlight is, from (A16) 

(Q/A) , = 0.0328 W/cm 2 or 104 BTU/hr ft 2 , 

rad a 

The calculated heat flux for the experimental conditions are 
summarized in the following table. 

The comparison between these, calculated heat fluxes and experi- 
mental data of Richards et al. >[1962] may be seen in Figure Al. 


Calculated steady state heat flux for the experimental 
conditions of Richards et al. [1962] 
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SPECIFIC HUMIDITY.gr/lb 
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cooled test section. 



A5. 2 Comparison with the Data of Ruccia and Mohr [i 960 ] 

Ruccia and Mohr present data for natural convection, and 
5.33 m/s (10 knot) and 12.6 m/s (24.4 knot) wind velocities. 

For natural convection over vertical flat plates McAdams [1942] 
recommends the following form of correlation (assumed applicable to 
vertical cylinders): 

Nu = a(Gr Pr) b . 

For L = 6.71 m (22 ft) and AT = 90 K (162°F), and the properties of 

1 2 

air of sec A2. 1 we get Gr Pr = 5. 22 (10) , which is beyond the range 

of the equation given by McAdams. We have fit the data given by 
McAdams at high Gr Pr by the following equation 

Nu = 0.0325 (Gr Pr)°* 384 (A21) 

from which we obtain 

(Q/A) = K/L Nu AT = 0.0751 W/cm 2 (238 BTU/hr ft 2 ) 

c 

and from (A5) 

( Q/A) dif = K/L Nu 3.38 (10) 3 H = 2.82 H W/cm 2 (H [ = ] g/g) 
or 

(Q/A) ... = 1. 275 H BTU/hr ft 2 (H [=] grains (lb), 
dii 

In Ruccia' s and Mohr's experiments the humidity varied from 

0.006 g/g (42 grains /lb) to 0.014 g/g (98 grains/lb). 

The forced convection Reynolds numbers corresponding to the 

5 5 

10 and 24 knot points are respectively 1.98 (10) and 4.83 (10) . 

These fall within the range of equations (A14) and (A15). From these 
equations 
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2 0 982 

(Q/A)' BTU/hr ft = 42.4 V ’ (V[ = ]ft/s) (A22) 

or 

2 - 0982 

( Q/A) dif BTU/hr ft = 0.2272 V * H (V[ = ]ft/s, 

H [ = ] grains /lb). (A 23) 

The radiation heat flux is, from (Al6), assuming T = 540°R 

(Q / A ) . = 101 BTU/hr ft 2 

Rad 

2 _4 0 982 

(Q / A) w / cm =4. 25 (10) V ’ (V [ = ] cm/s) 
c 

(Q/A) W/cm 2 = 0. 123 V 0,982 H (V [=] cm/s, H = g/g). 


The calculations are summarized in Table A2 and compared with 
the experimental data in Figure A2. 

Table A2 


V V H H 

cm/s ft/s g/g grains/lb 
Natural Convection 

(Q/A) 

c 

(Q/A) (Q/A)_ 

dif Rad 

W/cm 2 

< Q/i W 

(Q/A, total 
BTU/hr ft 

- 0. 006 

42 

0. 0751 

. 0170 

0. 0318 

0. 124 

393 

- 0.014 . 

98 

0. 0751 

. 0394 


0. 146 

464 

Forced Convection 







533 17. 5 0. 006 

42 

0. 223 

. 0501 


0. 305 

967 

1260 41. 2 0. 006 

42 

0. 515 

. 116 


0. 663 

2101 

533 17. 5 0. 014 

98 

0. 223 

. 117 


0. 372 

1179 

1260 41.2 0. 014 

98 

0. 515 

. 270 


0. 817 

2591 
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HEAT FLUX, W/cm 


WIND VELOCITY, ft/i 



Figure A2. Comparison of calculated steady state heat flux with the 
experimental data of Ruccia and Mohr [I960] for a liquid oxygen cooled 
vertical cylinder 45. 7 cm (18 in) OD x 6. 71 m (22 ft) long. 
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HEAT FLUX, BTU/h ft 




A5 . 5 Comparison with Data of Holten [1961] 

For natural convection heat transfer to spheres Bird et al. [1964] 
recommends 

Nu = 2 + 0.6 (Gr) 1/4 (Pr) 1/3 . (A24) 

The constants are 

D = 13.6 cm 4- 0.6 cm (approximate frost thickness x 2) 

= 14. 2 cm (or 0. 469 ft). 

AT = 90 K (or 162° F). 

(Other constants from Sec. A2. 1). 

From (A24) and the constants, 

h = 8. 16 (10)~ 4 W/cm 2 .K{or 1.44 BTU/hr-ft 2 - °F) (A25) 

(Q/A) = 0.0735 W/cm 2 (or 233 BTU/hr-ft 2 ). 

From (A25) and (A5), 

(Q/A) dif W/cm 2 = 3.38 (10) 3 hH = 2.75 H (H [ = ] g/g) 
or 

(Q/A) BTU/hr-ft 2 = 1. 25 H (H [ = ] grains/lb). 

The radiant heat flux is, from Sec A5. 2 

(Q/A) Rad = 0.0319 W/cm 2 (or 101 BTU/hr-ft 2 ). 
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Calculated steady state heat flux compared to the experimented, values 


are as 

follows: 





H 

g/g 

H 

grains /lb 

(Q/A) 

Exp (Q/A *Calc 
-BTU/hr ft 

(Q/ A ) E:cp <Q/A) Calc 

0.0065 

45.5 

3 25 

391 

0. 102 

0.123 

0.0070 

49 

315 

395 

0.099 

0.125 

0.0074 

52 

. 330 

399 

0.104 

0. 126 

0.0076 

53 

340 

400 

0. 107 

0. 126 

0.0097 

68 

336 

419 

0. 106 

0. 132 

0.0120 

84 

345 

439 

0. 109 

J 

0. 138 


It is apparent that the calculated heat flux is 20 to 30% higher 
than the experimental data. A possible explanation is that the heat 
flux in Holten's apparatus was impeded to some degree by the sheltered 
placement of the spherical vessels. 

Since this over-prediction is not apparent in the other compari- 
sons we assume that the predictive method is satisfactory. 1 
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Appendix B. Cooldown Surge 
Bl. Introduction 

Large amplitude pressure surges are often observed when a 
cryogenic liquid is forced into a warm transfer line. While water 
hammer pulses (sec. 4. 2) are associated with the impingement of 
liquid on a closing valve or other restrictions, these cooldown surges 
are associated with the opening of an inlet valve or the starting of a 
pump. Cooldown surges have been analyzed by means of computer 

r** 

solutions of the unsteady two phase flow equations (Steward et al. [1968]), 
the result of which being a stepwise calculation of all the variables in 
time and space coordinates. Not only are such detailed calculations 
too lengthy for this system analysis, but they produce unneeded infor- 
mation. In the present case only the peak amplitude which results 
from a particular mode of operation is required; therefore, a correla- 
tion of peak surge pressure as a function of the operating conditions 
is desirable. The following sections describe the development of the 
correlation used in the present system analysis and show some of the 
results in comparison with previous experimental data. 

The usual dimensional analysis approach to the development of 
an empirical correlation assumes that a proportionality exists between 
powers of dimensionless groups of parameters. It does not normally 
allow for the possibility of additive terms or other functional relation- 
ships such as trigonometric, exponential, etc. In the present approach 
both the dimensionless grouping of parameters and the forms of the 
relationships are derived from a model of the process. It should be 
understood that the derivations are not rigorous and leave empirical 
coefficients to be established from the experimental data. 
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During cooldown of pipelines with length to diameter ratios of 
450 and greater, observed stream and wall temperatures followed a 
"cold front" pattern. That is, at a given time temperatures upstream 
from a moving cold front were near the liquid saturation point; down- 
stream from that point the temperatures rose sharply to near the 
initial temperature. The limiting extreme, in which the temperature 
distribution in the line is reduced to one moving temperature step, 
from liquid to ambient temperature has been used with good results 
as a basis for cooldown time and pressure surge analyses (Steward 
[1965] [1968]) and will be used here. 

B2. Derivation 

As liquid flows into the warm part of the pipe and, after a time 
delay, becomes saturated, vapor begins to form in the liquid stream. 

The vapor carries liquid with it as it moves downstream and thus 
further increases the heat transfer surface area. As the rate of vapor 
formation increases, the pressure in the liquid-vapor transition region 
(or at the assumed liquid warm gas interface) rises due to the resistance 
to the flow of gas. Experimental data show that the measured pressure 
is in phase with the gas velocity. Thus, the surge process is unlike a 
frictionless spring -mass oscillating system in which the spring force 
is 90° out of phase with the velocity. Therefore, we conclude that 
the inertial effects (time derivitives of velocity) are small compared 
to the frictional effects. This conclusion is supported by detailed 
computations using the complete momentum, energy, and mass con- 
servation equations (Steward et al. [1968]) in which the magnitude of 
individual terms can be examined. Experimental measurements and 
computations show that the peak pressure occurs approximately at 
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the interface between liquid and vapor. The compressible gas in the 

warm part of the pipe expands and accelerates in the direction of 

flow to a maximum velocity at the discharge. The discharge velocity 

may be subsonic or sonic; however, for high peak pressure conditions 

which would be of concern, the discharge velocity would be sonic. For 

example in Fanno flow with f L /D (friction factor times gas length 

§ 

to diameter ratio) greater than 15 a ratio of peak pressure to discharge 

pressure, P / P greater than 5.5 would produce sonic (choked) 
p d 

discharge . 

For choked Fanno flow (constant area ideal gas flow with 
friction, without heat transfer) Shapiro [1953] gives 


ik z i~ M ' 2 + 111 In [lllltL 

D K 2 K i ( i + r-i M *) 

2 


where M^is the upstream Mach number and y is the ratio of specific 
heats. . For the purpose of correlation the logarithm term may be 
eliminated with no more than 20% loss of accuracy. Then, 


M,= 


/ 




(Bl) 


The mass flow rate corresponding to is 


w = M p AV where A is pipe area 
1 g s 


(B2) 


Substitution of ideal gas relations for density p and speed of sound V , 

% ® 

and equation (Bl) for M^yields, 
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where 


R is the gas constant. 


T 

t 


is gas temperature, and 


P is the peak pressure. 

P 

The total gas flow rate w may be expressed as the sum of an 
equilibrium flow rate w , 3111 excess vapor formation rate w ex * 


• • • 
w = w + w 

eq ex 


(B4) 


The equilibrium flow rate w e q is that which would exist in a 
steady cooldown without surging. In such a steady cooldown process 
most of the pressure drop P - P d (Pi is the upstream pressure) appears 
across the gas stream until the liquid front very nearly reaches the discharge 
end, i. e. , until 


cooldown is nearly complete. Then, for 


S 


<< I 

f* 





(B5) 



(B6) 
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Equation (Bl) is not strictly correct when the pressure P^is (as it is 
likely to be) below the value which would produce choking. The error 
in using equation (Bl) for subsonic flow is not great, however, because 
the upstream Mach number for a given mass flow rate and fL/D depends 
only slightly on the discharge mach number. The following example 
illustrates this point: 

Let fL/D = 15 and the discharge mach number = 0. 57. 

Fanno flow calculations for Y = 1. 4 gives an upstream mach number 

of 0. 194. When the discharge mach number is sonic (M^ = 1), then 

M^ = 0. 197 -- only 1. 6% higher. However, the corresponding pressure 

ratios are quite different, being P, /P, = 3 for the subsonic and P/P, = 

Id Id 

5. 5 for the sonic case. Thus, equation (Bl) which indicates that the 
upstream mach number is independent of the pressure ratio is approxi- 
mately correct. 

The excess vapor formation rate w is that which would result 

ex 

from liquid flows into the warm pipe faster than the rate which would result 
in steady cooldown. Such would be the case if the entering liquid were 
subcooled; then vaporization, which limits the inflow, would be delayed. The 
excess vapor formation rate is: 


W/ 

ex 


oc 


& h v + FAH S 


(B7) 


where L is the added heat transfer length, Ah is the heat of vapori- 

q v 

zation and F is some fraction of the enthalpy change of the gas, Ah , 

s 

between the saturation temperature and ambient temperature. 

The heat transfer length is assumed proportional to a reference 
liquid velocity and the saturation delay time t^, 

L OCU, t ,. (B8) 

q i sd 
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The delay time, i. e. , the time required for a mass of fluid 
entering the warm region to become saturated, is 



where Ah is the subcooling enthalpy, q is the heat flux per unit 
sc 

surface area andL^is a length corresponding to the excess mass of 
fluid entering, l ' is an effective heat transfer length during the time 

6 X 

t If we assume L / L ' - C ' then 

sd ex ex 



(B9) 


We assume the reference liquid velocity is proportional to 
the velocity of liquid before saturation, the vapor and gas in the pipe 
merely being displaced by the liquid. An additional restriction is 
imposed by a valve which may be in the process of opening. Alternately, 
a pump starting up may be thought of as a constant pressure source 
followed by a variable restriction. 



R v = o -t v i 

2 

K is a restriction factor. 

Substitution of (B8) and (B9) into (B7) gives 


Cj - c Ahis, p 
««' A tFAh, * 

where C is a new constant of proportionality which includes C' . 



(Bll) 
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As sume ; 


4- 


f, 


1 


c '(£)’ 


(B 12) 


where — is the pipe relative roughness 


1 


For an ideal gas p = — — 

g R r 

From (BIO), (Bll) and (B12) and (B13) we get 


<& 2 &) 

frP d )P, 


c, ( A + -yS) 


1 


[ B 1 3) 


(B 1 4) 



analyses for nitrogen and hydrogen as 0. 43. This value was considered 

constant. Also from the computer analysis liquid and gas friction factors 

e 

vary from 0. 01 to 0. 014 for an apparatus in which — ^ 0. 0001. There- 
fore, assuming an average f of 0. 012, we get C^, and = 0. 26. 

From liquid nitrogen experimental data for short valve opening 

time (t < t ,/2) 
v sd 

C = 0. 326 

From liquid nitrogen data with 16 second valve opening time 
2K = 0. 624 
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Comparisons of peak pressures from equation (B14) (solid, 
and dashed lines) with experimental data points are shown in figures 
B1 through B5. In all cases, the experimental points symbolized by 
circles are for various degrees of subcooling with highest peak 
pressures corresponding to maximum subcooling (storage temperatur 
at atmospheric boiling point). The solid lines also correspond to 
maximum subcooling in equation (B14). 

The square symbols and dashed lines denote saturated inlet 
liquid. It may be seen from these figures that the correlation gives 
upper and lower limits of peak pressure which envelope eventually ali 
of the experimental data as desired. 



a 40 



Experimental data, subcooled 
and partially saturated liquid 
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□ Experimental data, 
saturated liquid 

Equation B (14), 
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INLET PRESSURE RATIO, P/Pj 


Figure Bl. Peak surge pressure for liquid nitrogen in a 61 m (200 ft), 
1. 6 cm (0. 625 in) OD copper tube transfer line. 
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INLET PRESSURE RATIO, P,/P d 


Figure B2. Peak surge pressure for subcooled nitrogen in a 46 m (150 ft), 
1. 6 cm (0. 625 in) OD copper transfer line. 
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o 



INLET PRESSURE RATIO, P,/P d 

Figure B3. Peak surge pressure for liquid nitrogen in a 25 m (82 ft), 
1. 6 cm (0. 625 in) OD copper transfer line. 
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PEAK SU 


O Experimental data, subcooled 


and partially saturated liquid. 
□ Experimental data, 



INLET PRESSURE RATIO, P|/Pd 


Figure B5. Peak surge pressure for liquid hydrogen in a 61 m (200 ft), 
I. 6 cm (0. 625 in) OD copper transfer line. 
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Appendix C. Loading of the Minicomputer 


The method of loading programs into the minicomputer is described 
in Section 3.3 of the manual "How to Use the NOVA Computers" (1971). 
This appendix is an elaboration of the procedure described there. 

The Binary Loader 

The binary loader is usually in the minicomputer and is not 
destroyed by any of the programs. For that reason, following the 
procedure below should be sufficient to load any program in the mini- 
computer. 

However, occasionally a sharp spike on the power lines or a 
malfunction of the memory in the computer can cause the binary loader 
to be destroyed. If this should happen, the "Bootstrapping" procedure 
must be performed first. 

1 . Preparing the teletype 

The minicomputer and CRT display must be turned on but not 
running and the CRT display LOCAL/ LINE switch set to LOCAL. The 
teletype is plugged into the back of the CRT display. Turn the teletype 
power switch, LINE/OFF/LOCAL, to LINE. With the teletype reader 
switch in OFF position, place the paper tape containing the program 
in the reader with the tape coming out toward the front of the teletype 
and the narrower side from the sprocket holes in the tape toward the 
outside of the teletype. Position the tape so that unpunched tape is over 
the read position. Set the reader switch to ON. 
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2 . Loading the Program 

On the minicomputer, the upper row of switches is called the 
Switch Register, and are labeled 0 through 15. The lower row of 
switches are all momentary switches used for control and their 
functions are listed on the panel. 

Set the Switch Register 0, 1, 2 down and 3 thru 15 up. Press 
control switch RESET and then START. The teletype should start 
reading paper tape. At the end of each binary block on the tape, the 
printer mechanism should receive some null character and, 
thus, cycle or "chuckle 1 1 . Should it fail to do this, back up to the 
beginning of the paper tape and try again. Should it still fail it 
may be necessary to follow the "Bootstrapping 11 procedure 
and then try again. If the program loads correctly, it may self 
, start or it may halt the computer. If the computer halts, set 
all the Switch Register switches down except 14, which must be up, 
and press START. If subsequent operations are to be between the 
minicomputer and the CRT display, the LOCAL /LINE switch on 
the CRT should be set to LINE before pressing START, and the 
teletype turned to OFF. 

If the computer halts in the middle of the tape reading, the 
teletype has made an error while reading the tape. It is not necessary 
to start the tape over completely. Merely back the tape up at least 
one block and start the computer again. The beginning of a block 
is fairly easy to recognize. It consists of two frames with no punches 
in the tape except the sprocket hole, followed by a frame with at least 
the four top most positions punched followed by a frame with all 
positions punched. Reposition the tape on one of the two unpunched 
positions, set the teletype reader to ON and push START on the 
minicomputer. If the teletype fails to start "Chuckling", the tape 
should be restarted at the beginning. 
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Bootstrapping 

It may occasionally be necessary to restart the minicomputer 
from a completely obliterated condition. To do this, turn on the 
minicomputer and the CRT display. Set the CRT display LOCAL/ LINE 
switch to LOCAL. Prepare the teletype with the tape marked 
''CHUCKLING LOADER" in the reader just as you would for any 
other program tape. Next, follow the instructions in TABLE A-l. 
After the last line of the table has been executed the computer should 
read the tape and stop at the end. It will not "chuckle" the teletype. 

If it does not, read the tape the computer or teletype is likely in 
need of repair. If it does read the tape correctly, proceed to 
load program as described in the section. Binary Loader. 
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Set Switch Register Switch Position: ( U = up; D = down) Then Press 
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